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ABSTRACT 

We present Far Ultraviolet Spectroscopic Explorer (FUSE) and X-ray Multi- 
mirror Mission {XMM-Newton) data for the reverse shock of the 0-rich super- 
nova remnant (SNR) IE 0102.2-7219 in the Small Magellanic Cloud (SMC). The 
F?75'£' observations cover three regions with significantly different optical [Olll] 
intensities, all associated with the relatively bright part of the X-ray ring. Emis- 
sion lines of O VI AA1032, 1038 are clearly detected in the FUSE spectra. By 
combining this O VI doublet emission with the O vil triplet and O vill Lya fluxes 
from the XMM-Newton spectra and assuming a non-equilibrium ionization (NEI) 
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model with a single ionization timescale for the spectra, we are able to find a 
narrow range of temperatures and ionization timescales that are consistent with 
the respective line ratios. The increase of the ionization timescale r from North 
(r ^ 0.6 X 10-*^^ s cm"^) to South-East (r ^ 2 x 10-*^^ s cm~^) is indicative of 
increasing density in the X-ray bright ring, in good agreement with the optical 
[Olll] emission which is strongest in the South-East. However, if we assume a 
plane-parallel shock model with a distribution of ionization timescales, the O VI 
emission appears to be inconsistent with O Vll and O Vlll in X-rays. 

We also analyze the total XMM-Newton EPIC-MOS 1/2 spectra for the three 
regions. The X-ray spectra are dominated by strong emission fines of O, Ne, and 
Mg, however, wc detect an emission component that accounts for 14 - 25% of the 
flux and can be attributed to shocked ISM. We find that there is no consistent 
set of values for the temperature and ionization timescale which can explain 
the observed line ratios for O, Nc, and Mg. This would be consistent with a 
structured distribution of the ejecta as the O, Ne, Mg would have interacted 
with the reverse shock at different times. 

Subject headings: Shock waves — supernova remnants — ultraviolet: ISM — 
X-rays: individual (IE 0102.2-7219) 

1. Introduction 

The supernova remnant (SNR) IE 0102.2-7219 (hereafter E0102) is the brightest X-ray 
SNR in the Small Magellanic Cloud (SMC) and was discovered with Einstein (Seward & 
Mitchell 1981). In X-rays, it has a ring-like appearance and its emission is dominated by 
H-likc and He-like ions of O and Ne (e.g., Gactz ct al. 2000; Hughes ct al. 2000; Rasmussen 
et al. 2001; Sasaki et al. 2001; Flanagan et al. 2004). It consists of a bright X-ray ring and 
a much fainter plateau with a sharp outer edge that surrounds the ring (Fig. 1). This bright 
X-ray ring is caused by the reverse shock of the SNR, while the outer boundary of the faint 
plateau marks the location of the blast wave. The diameter of the bright ring is ~28", and 
the diameter of the outer rim is ~44". The radio emission in the 6 cm wavelength band 
(Amy & Ball 1993) is located outside the bright X-ray emission but within the outer faint 
X-ray rim (Gaetz et al. 2000). 

In the optical band, there is diffuse Ha emission surrounding the remnant, but a dis- 
tinctive Ha 'hole' at the location of the SNR. A filamentary shell of optical [O ill] emission 
with a diameter of ~30" is seen (Dopita et al. 1981), which is presumably emission from 
the ejecta. Further outwards, the remnant is surrounded by a faint, diffuse, partial shell of 
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[O III] with a diameter of ~60" that is thought to be the result of photoionization by the 
precursor. Emission Uncs of O, Ne, C, and Mg have been detected in the UV and optical 
spectra (Dopita et al. 1981; Morse et al. 1998; Blair et al. 2000), but H, He, and Fe are no- 
tably absent. Tuohy & Dopita (1983) estimated an age of ~1000 yrs for the remnant based 
on the optical diameter and the line of sight velocity range. Based on abundance estimates 
from optical and UV data, Blair et al. (2000) suggested that the progenitor was possibly a 
massive O star that, after becoming a Wolf-Rayet star, exploded as a Type lb supernova. 
More recently. Chevalier (2005) argued that E0102 is the remnant of a Type IIL/b supernova 
with a progenitor that lost most of the H envelope before the explosion. 

We present new observations with the Far Ultraviolet Spectroscopic Explorer {FUSE) 
to study the Ovi emission in regions in the X-ray bright ring with differing distributions 
of [Olll] and X-rays. In addition, we analyze X-ray observations with XMM-Newton and 
Chandra and compare the O emission in the Far-UV (FUV) to the X-ray emission. 

2. Observations and Reductions 
2.1. FUSE Data 

The FUSE satellite was launched in June 1999; details of its operation and on-orbit 

performance arc provided in the papers by Moos et al. (2000) and Sahnow et al. (2000). 
The FUSE instruments cover the 905 - 1187 A band with a nominal point source spectral 
resolution of i? ~ 20, 000. The FUSE telescope has four optical channels, each with its own 
mirror and grating. Two mirrors and gratings are coated with silicon carbide (SiC channels), 
and the other two mirrors and gratings are coated with Al and a lithium fluoride overcoat 
(LiF channels). One SiC/LiF channel pair is sensed by one microchannel plate intensified 
detector, and the other channel pair is sensed by a second detector. The SiC channels provide 
wavelength coverage from ~905 A to 1100 A, whereas the LiF channels cover the wavelength 
range of ~1000 A to 1187 A. Each of the two detectors is divided into two segments (lA 
and IB, 2A and 2B), so that 8 spectra are obtained for each pointing (LiFlA/B, SiClA/B, 
LiF2A/B, SiC2A/B). 

At the focal plane of each channel is an assembly holding three spectrograph apertures, 
dubbed LWRS, MDRS, and HIRS. In time-tag data mode, spectra through all three apertures 
are recorded simultaneously, but only one aperture is designated as primary. The primary 
science apertures for our program are the 4" x 20" MDRS apertures. 

After launch, small thermally-induced distortions were discovered that cause misalign- 
ments of these apertures. While periodic alignment activities are performed to keep the 
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channels approximately aligned, motions of 3 - 8" occur on an orbital timescale (~100 min- 
utes), preventing detailed co-alignment during a given integration. The LiFl channel was 
used for guiding in the data presented here, so that the location of the source is fixed for 
this channel. Due to the drift, not all the counts from the source fall on the detector. As 
this is not taken into account when the flux of the calibrated spectrum is computed, the 
flux of the SiC and LiF2 channels become lower than the real flux which is measured in the 
LiFl channel. The spectra of the SiC and LiF2 channels can be corrected by scaling the 
continuum flux to the correct continuum flux as determined in the LiFl spectrum. 

We observed three positions around the bright X-ray ring of E0102. The positions were 
selected to have differing distributions of [O iii] in the optical and O vii, O viii in the X-ray 
(Table 1). The observations were scheduled at times where the aperture position angles 
could be tailored to the positions being observed. Figure 1 shows the projected positions on 
a Chandra image and on an HST MVFPC2 [Olll] image. The net FUSE integration times 
are 12.4 ksec for the north pointing (N), 15.6 ksec for the northeast pointing (NE), and 8.2 
-|- 5.2 ksec for the southeast pointing (SE, two visits). We also performed a dedicated off- 
source position with MDRS for estimating the background. However, the spectra from the 
background pointing show more signiflcant contamination by background starlight than the 
on-source positions. Therefore, the background observation has not been subtracted from 
the source spectra. 

The data have been processed by the CalFUSE 3.0 pipeline, and the resulting spectra 
are shown in Figure 2. The only emission lines detected from the SNR are the broadened 
O VI doublet lines. We have modeled these spectra using functions available in IDL and 
routines from an IDL library including the fltting function MPFIT^ that performs least- 
squares curve fltting. Line fluxes and derived surface brightnesses in Ovi are provided in 
Table 2 (see Sect. 3.1 for more details). 

2.2. X-ray Data 

There are numerous public archival X-ray data sets for E0102 obtained with Chandra 
(Weisskopf et al. 1996) and XMM-Newton (Jansen et al. 2001). We flrst considered using the 
Advanced CCD Imaging Spectrometer (ACIS) data from the Chandra archive. Chandra data 
have far superior spatial resolution in comparison with XMM-Newton. However, we have 
determined that the brightest portions of E0102 as observed with the ACIS-S array in full 
frame mode are affected by pile-up, while the ACIS-S sub-array data and ACIS-I array data 
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with the remnant on-axis have poor photon statistics. Pile-up causes two or more photons to 
be detected as a single event. As the observed energy is the sum of the energies of two or more 
photons, the energy spectrum is thus distorted. Compared to Chandra observations, XMM- 
Newton observations arc less affected by pile-up due to shorter frame-times and lower spatial 
resolution (thus smearing emission from bright sources over a larger region of the detector). 
XMM-Newton is equipped with three European Photon Imaging Cameras (EPIC); one is a 
pn-type CCD (EPIC-PN) and two are MOS-type CCDs (EPIC-MOSl/2). Both the ACIS 
and the EPICs are sensitive in the energy band of ~0.3 - 10.0 keV. For SNRs, we are in 
general interested in soft X-rays (< 4.0 keV). Because the sensitivity of XMM-Newton in the 
energy range of 0.5 - 2.0 keV is higher and the spectral resolution of EPIC-MOS detectors 
is better than that of ACIS. wc decided to use XMM-Newton data primarily, except for the 
examination of the faint outer rim. The outer rim is much fainter than the bright ring and 
has much less pile-up. Chandra^s spatial resolution helps to isolate the rim emission from 
the bright ring emission. 

We selected archival data of three XMM-Newton observations which were taken early in 
the mission (April 2000 and April 2001): ObslDs 0123110201, 0123110301, and 0135720601. 
Because we want to study the emission lines, we use EPIC-MOSl and M0S2 data which have 
better spectral resolution than EPlC-PN. The net exposure times after removing times with 
soft proton flares are ~16 ks, 13 ks, and 20 ks for the observations 0123110201, 0123110301, 
0135720601, respectively. The data were obtained in small window mode resulting in a short 
frame time (0.3 s; for comparison, frame time of the Chandra ACIS full frame mode is 3.2 s) 
and negligible pile-up. The data are processed with XMMSAS 6.1.0. We extracted spectra 
of the reverse shock at the positions of the FUSE apertures (Fig. 3). The brightest lines in 
the X-ray spectra at these positions are the O Vll triplet lines and the O Vlll Lya line. 



3. Analysis 
3.1. O VI Doublet 

Although the diffuse O vi emission lines are broad and faint, these line features are well 
detected in the FUSE LiFlA and LiF2B spectra. The SiClA and SiC2B spectra also cover 
the spectral range around the O VI lines, however, the statistics of those data are poor due to 
the lower effective areas. Also, the SiC channels are more affected by the misalignment of the 
channels than the LiF2B channel and the emission from the brighter parts of the remnant 
may have largely moved out of the aperture. We see significant O VI emission in both LiFlA 
and LiF2B data. The rescahng of the LiF2B spectra by effective area and continuum flux 
shows that the LiF2B seems to have mainly covered the same O VI emission as the guiding 
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LiFlA channel. Therefore, we combine the LiFlA and LiF2B spectra for each pointing and 

analyze the co-added spectra around the O VI lines (1020 to 1050A). The strongest emission 
lines in the spectra of the E0102 pointings are terrestrial airglow lines. In the wavelength 
band around the O VI doublet AA1032, 1038 (1031. 926A, 1037.617A, respectively), there is a 
strong Hi airglow line at 1025. 7A, while the O VI emission is relatively faint, in part because 
of the broadening of the lines. We use this H I airglow line to calibrate the wavelength and 
to confirm that the spectral resolution is consistent with the aperture size. The O VI doublet 
fiuxes are derived by modehng the emission features between ~1030A and 1040A. 

Interstellar absorption is significant in the Far-UV and needs be taken into account 
when determining the line fiuxes. We include absorption lines of H2 and C 11 and correct 
the fiux for reddening through extinction. The H2 absorption has been derived using H2ools 
developed by McCandliss (2003). For the column density of H2 we use the Milky Way 
(MW) and SMC values from Tumhnson et al. (2002), i.e., median values from 13 stars in 
the northern part of the SMC. There is an evident absorption feature at ~1036A which is a 
Cll AA1036.3367, 1037.0182 doublet. We model this feature as a line at -1036 A with free 
width. We also model absorption by O VI along the line of sight, including both the MW 
and the SMC components. 

Wc combine night and day data (instead of using only night data to reduce the airglow 
emission) because of the poor photon statistics. Nevertheless, the sensitivity of FUSE allows 
the broad Ovi emission lines from this faint diffuse source to be well detected (see Fig. 2). 
The spectra arc binned with 0.3 A bins (corresponding to ~87 km s~^). Using the fitting 
routine MPFIT, we fit the Ovi AA1032, 1038 lines with two broad Gaussians, and the 
continuum is modeled as a constant. The widths of Ovi AA1032, 1038 Gaussian fines are 
tied together, and the relative positions of the lines are fixed to the wavelength separation: 
i.e. we fit the Ovi A1032 line (free wavelength) and tie the Ovi A1038 fine to the Ovi 
A1032 wavelength + 5.69lA. We also considered fitting each line with many zero-width 
Gaussians which fill the observed width of the broad lines, to see if we can find different 
velocity components in the lines. The zero- width Gaussians form pairs with fixed separation 
between the O VI A1032 and the Ovi A1038 components. However, the fits did not improve 
significantly. The normalization of each line is a free parameter except for the N pointing 
data, which have much poorer statistics. In the case of the N data, the normalization of the 
A1038 fine is fixed at half of the A1032 fine normafization (corresponding to the optically 
thin limit). From the line energy and width, we derive a velocity shift and FWHM. The 
fiux is corrected for extinction, using the Galactic extinction curve of Fitzpatrick & Massa 
(1990) and Fitzpatrick (1999) as well as extinction curves of the Magellanic Clouds from 
Gordon ct al. (2003). To obtain the de-reddened fiux, we use E{B — V) = 0.08, assuming 
that the extinction is caused half by the SMC and half by the MW, as suggested by Blair 
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et al. (1989, 2000). We compute the 90% confidence errors for the hne energy, width, and 
flux. The FWHM, shifts, and fluxes arc listed in Table 2. Figure 2 shows the spectra around 
A1036 with the fits overplotted and fit residuals. 

The Ovi emission is faintest in the N pointing and brightest in the SE pointing. In 
both the NE and SE pointings, the fits show that the gas has a velocity dispersion of about 

900 km s~^, while the faint Ovi line feature in the N pointing spectrum has a FWHM of 
1540 ± 350 km s^^. This is most likely not thermal line broadening, but instead indicates 
a velocity structure in the remnant (see Sect. 4.1). The rest SMC velocity is ~165 km 
(Danforth et al. 2002), so the broad Ovi emission lines are blue-shifted by ~100 km s~^ 
in the NE pointing and by ~330 km s^^ in the SE pointing. However, in the N pointing, 
the broad Ovi emission line feature is red-shifted by ~220 km s~^. This might indicate a 
non-symmetric expansion of the remnant. 

3.2. X-ray Emission Lines 

To combine the X-ray data with the FUV data, we perform an X-ray emission line 
analysis using the H-like and He-like lines of O, Ne, and Mg. The X-ray spectra have been 
binned with a minimum of 20 counts per bin to allow the use of statistics. We fit the 
EPIC-MOSl/2 spectra using the X-ray spectral fitting package XSPEC. The model consists 
of Gaussian hues at the known line energies of O, Ne, Mg, and Si (i.e. Ovil triplet at 
~0.57 keV, Ovill Lya at 0.654 keV, Neix triplet at ~0.91 keV, Nex Lja at 1.022 keV, 
Mgxi triplet at ~1.34 keV, Mgxil Lya at 1.473 keV, Sixill triplet at --1.84 keV, Sixiv 
Lya at 2.006 keV), a continuum-only model using the Astrophysical Plasma Emission Code 
(APEC, Smith et al. 2001), and line-of-sight absorption. The foreground absorption model 
consists of a component for the foreground MW absorption with a fixed A'h of 5.36 x 10^° cm~^ 
(Dickey & Lockman 1990) and a component with SMC abundances (Russell & Dopita 1992) 
to be fit. All the Lya hues are well fit with zero-width Gaussians. The CCD resolution does 
not allow the components of the triplet hues (O Vll, Neix, and Mgxi) to be resolved, but at 
least two zero-width Gaussians are necessary to reproduce the hne features. Therefore, we 
fit the triplet emission with two Gaussians at the energies of the forbidden and the resonance 
lines which are the two stronger lines and span the energy range of the triplet. The forbidden 
and the resonance lines are fit separately, i.e., their flux ratio is free in the fits, and the flux 
of the intcrcombination line is included in the fluxes of the fltted lines. We use the sum of 
the fluxes as the flux of the triplet. 

Oxygen lines are the strongest lines in both the FUSE and the XMM-Newton spectra of 
E0102. We compare line fluxes obtained from FUSE (Ovi doublet) and EPIC-MOS (Ovil 
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triplet and Ovill Lya) to line fluxes predicted by models. The observed line ratios are 
listed in Table 3. For the comparison with the observed lines, we use line fluxes computed 
with neiline^, which evaluates line emissivities for an ionizing or recombining plasma. 
Line fluxes are calculated by evolving the plasma at a given initial temperature for a given 
ionization timescale. As this calculation assumes a single ionization timescale for the plasma, 
we also model line fluxes using a plane-parallel shock model with a single temperature and 
a distribution of ionization timescales appropriate for a plane-parallel shock (VPSHOCK 
model in XSPEC; see Borkowski et al. 2001). We assume SMC abundances for all elements 
except for O, Ne, and Mg for both models. For the elements O, Ne, and Mg with enhanced 
abundances we use the numbers obtained from the global flts of the X-ray spectrum (see 
Sect. 3.3). 

The diagrams in Figure 4 show the hue ratios Ovil triplet/0 Vlll Lya, Ovi dou- 
blet/0 Vll triplet, and O VI doublet/0 Vlll Lya derived from FUSE scad XMM-Newton spec- 
tra and the error estimates based on the 90% confldence intervals for the line fluxes as a 
function of the ionization timescale r = Uet and temperature kT for the two models. In 
the VPSHOCK model, the parameter r is the upper limit of the ionization timescale that 
has a distribution characteristic for a plasma shocked by a steady plane-parallel shock. For 
the three line ratios for O, we evaluate confidence regions for the ratios as a function of kT 
and r. The lines in black show the ratio between O Vll and O Vlll observed in X-rays (O Vll 
triplet / O VIII Lya) . The cyan and orange lines compare O VI in Far-UV with O Vll and O Vlll 
( O VI doublet/0 VII triplet and O VI doublet/0 Vlll Lya). The results show that in the case 
of the neiline model, the confldence regions for the ratios Ovil triplet/0 Vlll Lya, Ovi 
doublet /O VII triplet, and Ovi doublet /O Vlll Lya intersect, yielding a well-defined kT vs. 
r = Uet region for each aperture: for N, kT ^ 0.3 — 0.4 keV and r ^ (0.5 — 2) x 10^^ s cm^^, 
for NE, kT ^ 0.3 - 0.4 keV and r > 0.4 x lO^^ s cm^^ and for SE, kT ^ 0.25 - 0.35 keV 
and r > 0.7 x 10"^^ s cm^^. While the kT-r regions for NE and SE overlap, the result for N 
indicates higher kT and lower r. However, the diagrams for the plane-parallel shock model 
show no overlap between the confldence regions of the line ratios. 

We also computed and compared the line ratios of OviI triplet/0 VIII Lya, NelX 
triplet /Nex Lya, and Mgxi triplet/Mgxil Lya to check whether different elements might 
indicate different plasma conditions. The results are shown in Figure 5. As the diagrams 
show, for neiline, the confldence regions for O, Ne, and Mg have little overlap and there 
are no conditions where all three intersect. There seems to be no single set of temperatures 
and ionization timescales that simultaneously describes the O, Ne, and Mg ratios. This is 
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consistent with previous results. In analyzing ASCA data for E0102, Hayashi et al. (1994) 
have found it necessary to introduce separate plasma components for each element. Analysis 
of Chandra High Energy Transmission Grating Spectrometer (HETGS) data has shown that 
the dispersed line emission images have different radii for lines of different elements, suggest- 
ing progressive ionization by the reverse shock propagating through the ejecta (Flanagan 
et al. 2004). Another explanation for the lack of overlap between the ions is that the X-ray 
emission arises from a distribution of shocks with different velocities as suggested for SNR 
1987A by Zhekov et al. (2005). In the case of a plane-parallel shock model (VPSHOCK) the 
three confidence regions overlap for kT > 1 keV and r < 10^^ s cm~^. However, when 
the O VI hne results from FUSE are included, there is no region of common temperatures 
and ionization timescales which satisfies the O line ratios. This might be indicating that the 
bulk of the O VI emission arises from a different location with different plasma conditions 
than the location from which the bulk of the O Vll and O Vlll emission arises. 



3.3. Global Fitting of the X-ray Spectrum 

To analyze the global X-ray spectra of the three regions, we first fit the spectra in 
XSPEC using the VNEI model with NEI version 2.0 for a plasma in non-equilibrium ioniza- 
tion (NEI) with variable abundances. This model assumes a single ionization timescale for 
the emitting plasma and is a simplification compared to the plane-parallel shock model (see 
Borkowski et al. 2001). The NEI version 2.0 models use a newer atomic database (Astro- 
physical Plasma Emission Database, APED, Smith et al. 2001b) than the versions 1.0 and 
1.1 and reproduces more accurate line emissivities. However, they do not contain inner-shell 
processes. For the absorption we include two components: one component represents the 
foreground MW absorption with a fixed iVn of 5.36 x 10^*^ cm~^ and the other represents 
the SMC absorption. For each region, the EPIC-MOSl and M0S2 spectra for the three 
observations are fit simultaneously with a single model. 

Although we fit the spectra with various sets of initial parameters, models with a single 
VNEI component fail to fit the spectra adequately. In all cases, the model predicts a lower 
flux than the observed spectrum next to the Mgxi triplet line feature, exactly where the 
Mg XII Lja line would be expected. Because the fits for the X-ray bright ring are dominated 
by the strong O VII, O VIII, NelX, and NeX lines, we conclude that the elements O, Ne, and 
Mg might not be ionized to the same extent, with Mg being more highly excited than the 
other elements. 

We therefore try a model with additional VNEI components to represent different plasma 
components. First, we include a component for the ISM shocked by the blast wave: To 



-10- 



estimate this component, we use Chandra ACIS-S3 data (ObsIDs 5123 and 5130) and extract 
a spectrum in the X-ray plateau covering the outer blast wave, i.e. a ring centered on RA 
= 01^ 04™ 02.0" , Dec = -72° 01' 53.2" with inner radius of 19" and outer radius of 23" . 
This emission is very faint and not affected by pile-up. In the XMM-Newton EPIC-MOSl/2 
data it is not possible to separate clearly the emission from the blast wave and the ejecta 
because of the poor spatial resolution of XMM-Newton. Fitting the blast wave spectrum 
extracted from Chandra data with one VNEI component (absorption component as above) , 
we obtain: kT = 1.1 keV, r = 5.7 x 10^° s cm~^. Therefore, we fit the XMM-Newton spectra 
of parts of the bright X-ray ring with a model including a VNEI component for the ISM 
with kT and r fixed, and all the abundances set to the SMC ISM abundances (Russell & 
Dopita 1992). For the ejecta emission, we include two VNEI components: a component for 
the emission dominated by O, fixing the parameters for kT and r to the results from the 
line diagnostics (see Fig. 4) with all other element abundances set to zero, and a component 
accounting for the rest of the ejecta emission, consisting of Ne and Mg (with abundances 
for the other elements set to zero), with free kT and r. The same set of absorbing columns 
(Galactic and SMC as described above) is used for both components. The fits are better 
than with a single VNEI model, and the fit parameters are hsted in Table 4. The spectra 
of the ObsID 0123110201 data and the best fit model, obtained by a simultaneous fit of all 
spectra of the three ObsIDs, are shown in Figure 3. 



4. Discussion 

4.1. Results from the Spectral Analysis 

If the line broadening in the FUSE spectra were due to thermal motion, the FWHM 
of about 900 km s^^ would correspond to an unreasonably high temperature of 3x10^ K. 
The analysis of the EPIC-MOS data shows that the plasma temperature is about 10'' K. 
There are at least three scenarios that can explain this discrepancy. First, it is possible that 
the collisionless heating of the electrons behind the reverse shock front is negligible. Hughes 
et al. (2000) argue that this is the case at the outer blast wave and that significant particle 
acceleration is required there. In such a case, the oxygen ions, and thus also Ovi, have 
a much high temperature than the electrons which are responsible for the X-ray emission. 
Second, the electron temperature might have been underestimated by applying a single 
ionization timescale model (VNEI). Third, the broadening might be caused by turbulence or 
by bulk motions rather than by thermal motions. Flanagan et al. (2004) derive expansion 
velocities of about ±1000 km s~^ from Chandra HETGS data for the bright X-ray ring. 
The velocities measured from the HETGS data could be explained by geometric effects, if 
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E0102 has, e.g., a tilted, expanding ring geometry or a distorted shell with density variations. 
Flanagan et al. (2004) suggest that the difference in velocity shifts in different parts of the 
remnant might indicate that the southeastern part of the remnant shell is approaching and 
the northwestern part is receding. This is consistent with the velocity shifts obtained for the 
broad Ovi emission lines which show blue-shifts in the NE and SE pointings, whereas the 
lines in the N pointing indicate a red-shift (Sect. 3.1). 

The X-ray emission of the bright X-ray ring is mainly ejecta emission. Because the 
spectra in all three regions are not well fit by a single component NEl model, we have fit the 
spectra with two NEl model components with variable abundances of O (first component) 
and Ne and Mg (second component) for the emission from the ejecta shocked by the reverse 
shock, and an additional NET model to account for the emission from the ISM shocked by 
the blast wave. We calculate the unabsorbed flux (0.3 - 10.0 keV) of the ejecta and shocked 
ISM for each region and compare the flux fraction: The ratio between the surface brightness 
of the total ejecta emission and the ISM emission, S'ejecta/5'isM, is 6, 3, and 3 for regions 
N, NE, and SE, respectively. The surface brightness of the ISM component is S'ism ~ 5 x 
-|^g-i5 gj,g g-i arcsec"^ in all three regions. For comparison, the surface brightness of the 
emission of the shocked ISM in the outer rim is S'outemm = 4.7 x 10^^^ erg cm~^ s~^ arcsec"^, 
in good agreement with the estimated ISM component for the emission of the X-ray bright 
ring at the positions of the FUSE pointings. The differences between the total flux of the 
three regions in the X-ray bright ring seem to be caused by differences in the emission from 
the ejecta shocked by the reverse shock as well as by variations in absorption as obtained 
from the spectral fits. 

The line ratio diagrams in Figure 4a show that the Ovi seen in FUV and the Ovil 
and O Vll X-ray emission are compatible with the assumption of an ionizing plasma with a 
single ionization timescale r, as calculated by neiline. In this case, the line diagnostics of 
the O emission obtained with FUSE and XMM-Newton show that the ionization timescale 
of the O gas differs in the three regions, increasing from N to SE: r ~ 6 x 10^° s cm~^ in 
the N region, 1 x 10^^ s cm~^ in the NE region, and 2 x 10^^ s cm~^ in the SE region. The 
temperatures in the three regions are comparable: kT = 0.3 - 0.5 kcV. This result may 
indicate that the density in the ejecta is higher in the southeastern part of the X-ray bright 
ring, consistent with the optical data which show that the [O ill] emission is strongest in the 
SE. If we instead assume a plane-parallel shock with a distribution of ionization timescales 
using the VPSHOCK model, we do not obtain a consistent set of kT and the upper limit 
of T for the O VI, O VII, and O VIII fines. However, the diagrams in Figure 5b indicate that 
the flux ratios of emission lines of O, Ne, and Mg in X-ray, and thus also O Vll and O Vlll, 
are consistent with a common temperature and ionization timescale assuming the plane- 
parallel shock model. As the diagrams in Figure 4b compare O VI in FUV with the oxygen 
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emission in X-ray, the fact that O VI, O Vll, and O Vlll show no overlap might suggest that 
the oxygen hne emission in the FUV and the oxygen hne emission in X-rays together are not 
consistent with the assumed model. As copious [O ill] emission is observed in the optical, 
there must also be denser shocked material and the O VI emission might be related to it as 
well. Therefore, multiple shocked plasmas or a range of velocities are possibly involved in 
the FUV and X-ray emission of E0102. 

These results show that we find a few possible scenarios for the emission from this SNR. 
However, as we arc limited by the size of the aperture of FUSE and the spatial resolution of 
XMM-Newton wc arc not able to constrain whether the emission arises from a single shock 
or multiple shocks. Moreover, the shock can be ionizing or recombining and thus radiative. 
Although we have optical [O ill] images and X-ray images with O Vll and O Vlll emission of 
high spatial resolution, we lack spatial information on O VI. If it were possible to obtain O VI 
imaging data we would be able to compare the distribution of O in the remant and thus to 
determine if the shock is ionizing or recombining. Furthermore, hydrodynamical models for 
reverse shocks are necessary to understand the distribution of elements as well as physical 
parameters like kT or r in relatively young SNRs like E0102. 



4.2. Mass Estimates 

To derive the ejecta mass, we numerically estimate the volume of the SNR included 
in the slit, assuming that the SNR emission is from a thick shell. The outer radius of the 
shell is assumed to coincide with the outer edge of the bright X-ray ring: Rg = 16'.'0±1'.'0 
= 4.7±0.3 pc (for D = 60 kpc, van den Bergh 1999); the inner radius of the bright X-ray 
ring is taken to be rg = 12'.'0±1'.'0 = 3.5±0.3 pc. The resulting volume, V, of the observed 
emitting gas is the intersection between this thick shell and a projected rectangle with the 
slit size {A = 4"x20" = 6.77 pc^ = 6.45 x 10=^^ cm^). The results are hsted in Table 5. We 
assume an error of ~20% for the estimated volumes, as there are systematic uncertainties 
including the actual shape of the shell. 

The Chandra HETGS spectra show that the plasma mainly consists of O, Ne, and Mg 
(Flanagan et al. 2004). The EPIC-MOS fits also show that the ejecta emission results from 
O, Ne, and Mg. Therefore, we assume a metal- rich plasma consisting of O, Ne, and Mg 
excited to H-like, He-hke, or fully ionized states. 

The XSPEC VNEI emissivity normalization (Table 4) is K = 10-^^/{47rD^) J neUndV, 
where D is the distance to the source [cm], nn is the hydrogen number density [cm~^], and 
Ue is the electron number density [cm~^] . In the case of a plasma with high overabundance of 
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elements other than hydrogen, nn needs to be replaced by a reference density, ^ref, that is a 
sum over the number densities of the included elements, thereby taking the abundances into 
account. The electron density is ne,i = a\ni, where n-i is the density of element i. The factor 
cti depends on the assumed mean charge state for each clement; we assume a = 7, 9, 11, for 
O, Ne, and Mg, respectively. Assuming D = 60 kpc= 1.85 x 10^^ cm and uniform density 
in the thick shell, we estimate no, nNe, and nug for the two components for each region 
(see Table 5). Using these number densities, we obtain the O, Ne, and Mg mass for the 
entire remnant. We include a filling factor / to account for the unknown volume fraction of 
the X-ray emitting plasma. The densities are proportional to /~°'^. Assuming a thick shell 
with an outer radius of Rs — 4.7±0.6 pc and an inner radius of Tg = 3.5±0.6 pc, we obtain 
Mo = 4 ± 2 f-^ Mq, Mnc = 1.6 ± 0.4 f-^ Mq, Mug = 0.2 ± 0.1 f-^ Mq. If wc assume 
a thick ring with V = 6.6 x 10^"^ cm~^ as estimated by Flanagan et al. (2004), wc derive: 
Mo = 4 ± 1 Mq, Mnc = 1.5 ± 0.3 f-^ Mq, Mug = 0.2 ± 0.1 f-^ Mq. We estimate a 
few Mq in O, indicating a massive progenitor. 

If we assume a volume fraction of unity, the oxygen mass Mq is consistent with the result 
from Chandra HETGS by Flanagan et al. (2004) who also assumed / = 1. By comparing 
their Mo estimate with the models of Nomoto et al. (1997), Flanagan et al. (2004) suggest 
that the mass of the progenitor of E0102 was ~32 Mq. Blair et al. (2000) compare abundance 
ratios based on optical and UV lines to values derived by Nomoto et al. (1997) and conclude 
that the E0102 progenitor is best described by a 25 Mq model. Our mass estimates for O 
and Mg as well as the mass ratios agree well with the 25 - 30 Mq model (Nomoto et al. 
1997), while the Ne mass estimate is shghtly higher than the model prediction. Even for a 
filling factor as low as / 0.1, the mass estimates are consistent with a massive progenitor 
(Nomoto et al. 1997). 

In these calculations, there are still some uncertainties that we have to be aware of. As 
the reverse shock has not yet reached the center of the remnant, there is still unshocked 
material interior to the reverse shock. If the unshocked material has the same abundance 
distribution as the shocked material (i.e., well-mixed ejecta) the ratios will not be much al- 
tered, but the estimated masses would be lower limits. Furthermore, the ionization potential 
increases with the atomic number; the shock might be strong enough to ionize O all the way 
up to the He-like and H-like stages, but not as efficient in ionizing Mg. Therefore, a part of 
Mg could be in lower ionization stages and not emitting in X-rays. Finally, assuming a too 
high volume filling factor can result in overestimating the total masses. Numerical calcula- 
tions have shown that the ejecta shocked by the reverse shock form thin layers that can be 
deformed by Rayleigh- Taylor instabilities, decreasing the fraction of the observed ejecta in 
the apertures (Chevalier et al. 1992; Blondin & Ellison 2001). 
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5. Summciry 

We have observed three regions in the bright X-ray ring of the 0-rich SNR IE 0102. 2- 
7219 in the SMC using the FUSE MDRS (4"x20") aperture. These three regions cover 
parts of the remnant ring with varying distributions of [O III] in the optical, and O VII and 
Ovill in the X-ray. The ring is caused by the reverse shock propagating into the ejecta. In 
the northern region (N), the FUSE spectrum shows very faint but velocity broadened O VI 
emission. In the northeastern (NE) and the southeastern (SE) parts of the SNR ring, there is 
significant O VI emission with broad O VI AA1032, 1038 lines, indicating a velocity dispersion 
of about 900 km s^^. The emission from the SE pointing appears to be blue-shifted by 
~200 km s~^ relative to the emission from the NE pointing. 

We have also analyzed the XMM-Newton EPIC-MOSl/2 spectra of E0102 in the same 
regions as in the FUSE observations. For the combined analysis of the FUV and X-ray 
data, we consider two cases: 1. shocked plasma described by a single temperature and 
single ionization timescale (using the models neiline and VNEI) and 2. plane-parallel shock 
with a distribution of ionization timescales (using the VPSHOCK model). In the first case, 
combining O VI doublet emission from FUSE data with O Vll triplet and O Vlll Lya in EPIC- 
MOSl/2 spectra yields kT ^ 0.4 keV and r 6 x 10^° s cm"^ for the N region, kT ^ 0.4 keV 
and T 1 X 10^^ s cm~^ for the NE region, and kT ^ 0.3 keV and r f« 2 x 10^^ s cm~^ 
for the SE region. The ionization timescales of the O gas seem to increase towards the 
south (tn < tne < tse), suggesting a higher density in the south. The total EPIC-MOSl/2 
spectra are well fit with a three component non-equilibrium ionization model; two VNEI 
components are necessary to describe the ejecta emission (one for the O plasma, one for 
Ne + Mg plasma) and an additional VNEI component is used to model the emission from 
the shocked ISM. The temperature of the Ne -|- Mg component is higher [kT > 1 keV) 
than the O component. For the plane-parallel shock model, however, the observed ratios of 
Ovi/Ovil and Ovi/Ovill are not consistent with the Ovil/Ovill ratio, in that there is 
not a common region of allowed temperatures and ionization timescales. When the triplet 
and Lya ratios are considered for O, Ne, and Mg in the X-ray bandpass alone, the observed 
ratios do have significant overlap in temperature {kT > 1.0 keV) and ionization timescale 
(r < 1 x 10"*^^ s cm~'^). These results might indicate that the FUSE and XMM-Newton data 
are sampling multiple shocks with different conditions. To better understand these results 
in the future, we still require information on the spatial distribution of O VI in FUV as well 
as hydrodynamic models describing the structure of the reverse shock. 

From the X-ray emission, we estimate the mass of O, Ne, and Mg in the entire remnant. 
Making the assumption that the remnant is a thick shell, we obtain Mq = 4 ± 2 f-^ Mq, 
= 1.6 ± 0.4 /°-^ Mq, Mug = 0.2 ± 0.1 /°-^ Mq. A comparison with previous measure- 
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merits and progenitor models shows that our results are consistent with a massive progenitor. 
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Fig. 1.— Chandra AC1S-S3 image (left) and HST WFPC2 [Olll] image (right) with FUSE 
MDRS aperture positions. See fl.jpg. 
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Fig. 2. — Close-ups of FUSE spectra of N, NE, and combined SE pointings. LiFlA and 
LiF2B spectra are combined. The thick dashed hne shows the model (see Sect. 3.1). The 
thick lines at the bottom mark the rest wavelengths of the O VI lines. The dotted line shows 
the H2 absorption in an arbitrary scale. 
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Fig. 3.— XMM-Newton EPIC-MOSl (solid) and M0S2 (dashed) spectra of the observation 
0123110201 for the regions N, NE, and SE. Fit model includes three VNEI spectral compo- 
nents (see Sect. 3.3). The prominent O, Ne, and Mg hnes are labeled in the N spectrum. 
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Fig. 4. — Line ratios Ovil triplet/0 Vlll Lya, Ovi doublet/0 Vll triplet, and Ovi dou- 
blet/0 vill Lya for the regions corresponding to the apertures of the FUSE pointings N, 
NE, and SE (solid line). The estimated range for the error is shown with dotted lines. To 
calculate the line ratios, we used the models a) neiline and b) VPSHOCK (in XSPEC). In 
the VPSHOCK model, r is the upper limit of the ionization timescale that has a distribution 
typical for a plane-parallel shock. 



-22- 



a) 



Reg 
EPIC 



ion N 
-MOS 



> 

(a 



plet/Lya 
'plet/lya 



Trif 
Ne Trif 

Mg Triplet'/Lya 




b) 



> 




Reg 

EPIC- 

Tnr 
Ne Tr , 

Mg Triplet'/Lya 



ion NE 

-MOS 

plet/Lya 
riplet/Lya 




10" 10 
T [s/cm^] 



Reg 
C 




ion SE 
MOS 



Trif 

Ne Triplet/Lya 
Mg Triplet/Lya 



plet/Lya 



> 





■\\y ■■, Region N 




\V\ / EPIC- MOS 




\ \VN\ ''O Triplet/Ly<3; - 
\ '■■■A--:-T\ Ne. Triplet/Lya 
\ ''■•\T\\ Mg'''-Triplet'/Lya 



T s/ cm 



Region NE 

EPIC-MOS 

Triplet/Lya 
Ne Triplet/Lya 
Mg Triplet/Lya, 




Region SE / 
EPIC-MOS--'" 

Triplet/Lya . 
Ne Triptet/Lya( 
Mg Tri/let/Lya" 



T [s/cm^] 




Fig. 5. — Line ratios Ovil triplet/0 Vlll Lya, Neix triplet/Nex Lya, and Mgxi 
triplet/MgXll Lya for the three regions corresponding to the apertures of the F?7>S'£' point- 
ings N, NE, and SE obtained from EPIC-MOS spectra (sohd hne, 90% errors as dotted 
hnes). The used models are a) neiline and b) VPSHOCK (in XSPEC). The parameter r 
of the VPSHOCK model is the upper limit of the ionization timescale. 
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Table 1. FUSE, XMM-Newton, and Chandra data used for the analysis. 



Obs. ID^ 


Pointing 


direction 


Aperture/ 


Start time (UT) 


Exposure 




RA 


Dec 


Instrument 




night -|- day'' 




(J2000.0) 






[ksec] 


FUSE 


C0830101 (N) 


01:04:01.74 


-72:01:39.1 


MORS 


2003-09-19 11:35:05 


12 


C0830201 (NE) 


01:04:04.24 


-72:01:45.1 


MDRS 


2002-06-06 23:51:41 


16 


C0830301 (SEl) 


01:04:04.04 


72:02:01.9 


MDRS 


2003-08-29 07:54:23 


8 


C0830302 (SE2) 


01:04:04.04 


-72:02:01.9 


MDRS 


2004-08-21 14:52:32 


5 


XMM-Newton 


0123110201 


01:03:47.25 


72:01:58.3 


EMl/2, SW= 


2000-04-16 20:09:26 


16 


0123110301 


01:03:47.19 


-72:01:58.3 


EMl/2, SW^ 


2000-04-17 04:43:36 


13 


0135720601 


01:03:45.79 


-72:00:51.1 


EMl/2, SW'^ 


2001-04-14 20:47:25 


20 


Chandra 


5123 


01:03:56.41 


-72:01:18.7 


ACIS-S 


2003-12-15 04:27:06 


20 


5130 


01:04:12.31 


-72:01:51.2 


ACIS-S 


2004-04-09 13:08:57 


19 



^Throughout the paper, the FUSE observations are called N, NE, and SE, as indicated in the table. 
^For FUSE data. 

'^EPIC-MOSl/2 data in small window mode. 
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Table 2. FUSE spectral fit results with one set of O VI AA1032, 1038 lines per each 

aperture position. 



Data 


Bin 


FWHM 


Line Shift 


Surface Brightness 


Reduced 










-Paio32 


-Pa1038 


x' 




[A] 


[kms-i] 


[kms-i] 


[ergcm~^ arcsec"^] 


[ergcm-^ s-i arcsec-^] 






0.3 


1540±350 


380±180 


(1.2±0.2) X 10-15 


(0.6±0.1) X 10-15 


1.59 


NE 


0.3 


970±120 


60±60 


(1.7±0.3) X 10-15 


(1.3±0.4) X 10-15 


1.65 


SE 


0.3 


830±310 


-160±30 


(2.3±1.1) X 10-15 


(2.0±1.1) X 10-15 


1.72 



Note. — The flux is de-reddened. The errors are 90% confidence ranges. 
'*A1038 line flux flxed to half of the A1032 flux. 
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Table 3. Line ratios for O, Ne, and Mg. 



Region 


VII triplet/ 


VI doublet/ 


VI doublet/ 


Neix triplet/ 


MgXI triplet/ 




VIII Lya 


VII triplet 


VIII Lya 


Nex Lya 


Mgxil Lya 


N 

NE 

SE 


1 4Q+0.93 


1.0412:1 


0.82l°:^« 
1 90+0.99 

-'-•^<J-0.90 

-1 -19+0.41 

-'-■-'-^-0.52 


1 40+0.26 

-■-■^^-0.20 

Z.UO_Q 24 

9 9c; +0.37 
^■■^ "J -0.20 


4.631?:^^ 

/ ■Do_6.96 



Note. — Emission line fluxes arc obtained from the FUSE spectra (O VI doublet) and 
from the XMM-Newton spectrum (He-like triplets and H-like Lya line of O, Ne, and Mg) by 
fitting Gaussians. 
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Table 4. XMM-Newton fit spectral results for regions corresponding to FUSE MDRS 

aperture positions. 





N NE 




SE 


A^H (SMC) [1021 cm-2] 


1.6(1.3-2.0) 0.1(0.0- 


0.4) 


0.4 (0.2 - 0.7) 


1. VNEI component (0 ejecta) 


fcTi [keV] 


0.40^ 0.38' 




0.30^ 


(solar^) 


93 (83 - 104) 86 (80 - 


94) 


99 (92 - 107) 


Ti [10^° s cm-3] 


6.0*^ 10.^ 




20.^ 


Ki' [10-6] 


2.4(2.2-2.6) 2.3(1.5- 


3.2) 


3.7 (3.5 - 3.9) 


2. VNEI component (Ne and Mg ejecta) 


m [keV] 


3.5(2.7-3.7) 1.5(1.0- 


3.2) 


1.1 (0.8 - 1.6) 


Ne (solar^) 


82 (75 - 89) 86 (66 - 


116) 


82 (77 - 94) 


Mg (solar*' ) 


24 (20 - 28) 33 (23 - 


49) 


30 (22 - 40) 


T2 [10^° s cm-3] 


1.8(1.7-1.9) 2.3(1.5- 


3.7) 


3.0 (2.1 - 4.4) 


K2' [10-6] 




= K, 


3. VNEI component (ISM) 


/v./ [10-^] 


3.0 (2.8 3.3) 2.0 (1.8 


2.. 3) 


2.3 (2.1 2.5) 


Fit information 




559.7 608.8 




618.8 


d.o.f 


374 368 




403 


Reduced 


1.50 1.65 




1.54 



Note. — The 90% confidence range is given in parentheses. The model consists of three 
VNEI components: one for the blast wave emission, two for the ejecta emission. The absorp- 
tion for SMC is modeled with SMC abundances. Additional Galactic absorption is modeled 
with fixed column density of 5.36 x 10^° cm-^. 

^Result from the line diagnostics (Fig. 4). 

''Abundances from Anders & Grevesse (1989). 

■^Normalization K = IO'^^/IAttD"^) J fienudV, where D is the distance to the source (cm), 
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Ue is the electron density (cm and nu is the hydrogen density (cm 
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Table 5. Estimated volumes observed through the MDRS apertures, number densities, 
and masses derived from the fits of EPIC-MO S spectra using a three component VNEI 

model. 





N 


NE 


SE 


V [pc^] 


19.8±4.0 


24.3±4.9 


26.0±5.2 


V [10'^^« em-^] 


5.8±1.2 


7.1±1.1 


7.6±1.5 


Number densities 


nop-^ [cm ^] 


0.045±0.009 


0.038±0.008 


0.050±0.009 




0.013±0.005 


0.012±0.004 


0.014±0.004 


umJ'-' [cm-3] 


0.0012±0.0004 


0.0014±0.0005 


0.0016±0.0005 


Mass^ assuming a thick shell 


Mo/f-' [Mo] 


4.3±1.6 


3.7±1.4 


4.8±1.7 




1.6±0.8 


1.4±0.7 


1.7±0.7 


Mujf-' [Mo] 


0.18±0.08 


0.21±0.09 


0.24±0.10 


Mass*^ assuming a thick ring*^ 


Mo/r-5 [Mo] 


4.0±0.8 


3.4±0.7 


4.4±0.8 


MNe//°-^ [Mo] 


1.5±0.5 


1.3±0.4 


1.6±0.5 


MMg//°-^ [Mo] 


0.16±0.06 


0.19±0.06 


0.22±0.06 



^O, Ne, and Mg masses in the entire remnant assuming the density in the regions N, NE, 
or SE. 

''Assuming a thick ring suggested by Flanagan et al. (2004). 
is the volume filling factor. 
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